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We show that the transition from Gaussian to the (/-Gaussian distributions occurring in atomic 
transport in dissipative optical lattices can be interpreted as self-organization by recourse to a 
modified version of Klimontovich's S-theorem. As a result, we find that self-organization is possible 
in the transition regime, only where the second moment (p 2 ) is finite. Therefore, the nonadditivity 
parameter q is confined within the range 1 < q < |, although whole spectrum of q values i.e., 
1 < q < 3, is considered theoretically possible. The range of q values obtained from the modified 
S-theorem is also confirmed by the experiments carried out by Douglas et al. [Phys. Rev. Lett. 96, 
110601 (2006)]. 



I. INTRODUCTION 



The importance of understanding physical phenomena exhibiting asymptotical inverse power law distributions 
becomes apparent if we consider its emergence in such diverse fields as subregion laser cooling [1], the heartbeat 
histograms of healthy patients [2], blinking nanocrystals [3], rheology of steady-state draining foams [4], econophysics 
[5], earthquake models [6], to name but a few. 

One such distribution is the Tsallis distribution [7] , which found many applications in different fields [8] . Recently it 
was pointed out that the momentum distribution of cold atoms in a dissipative optical lattice is a Tsallis distribution 
[9]. Soon, this theoretical prediction has been verified experimentally [10]. 

Three different regimes of atomic transport is possible in an optical lattice (see Ref. [11] for details), depending on 
the depth of the optical potential. In the intermediate regime between the diffusive motion in deep potentials and 
ballistic motion in shallow potentials, the atomic transport is governed by the Tsallis distribution [9, 10]. In fact, it 
has been shown, both theoretically and experimentally that a transition from an ordinary Gaussian to the g-Gaussian 
distribution occurs by changing the ratio of the recoil energy to the potential depth [9, 10]. 

Our aim in this paper is to show that this transition can be interpreted as a self-organization in the context 
of modified S-theorem. Klimontovich [12-17] originally proposed S-theorem (Self-organization theorem) in order to 
generalize Gibbs' theorem [18] to open systems, since the latter rests on the assumption that all the compared 
distributions have the same mean energy values. However, this assumption does not hold when one studies open 
systems where the mean energy of the system is not constant due to the interaction of the system with the environment 
(e.g., a laser field in the case of optical lattices). Therefore, Klimontovich defined a quantity called effective energy. 
Equating the mean effective energies of different states and then comparing the associated entropies, he showed that 
S-theorem orders the associated entropies in such a way that the state closer to equilibrium possesses a greater entropy 
compared to the other states. In other words, we have a more ordered state as the control parameter increases. This 
decrease of entropy on ordering is called self-organization by Haken [19]. We also mention that the process of equating 
mean effective energies of different states is called renormalization by Klimontovich. As a result, S-theorem did not 
only succeed in generalizing Gibbs' theorem but gained recognition as a criterion of self-organization in open systems, 
since it allows us to compare even the stationary nonequilibrium distributions [12-17]. 

We also mention that the use of S-thcorcm is not limited to analytical models. It has also been used for many 
numerical models such as logistic map [20], heart rate variability [21, 22] and the analysis of electroencephalograms 
of epilepsy patients [23] as a criterion of self-organization. 

However, S-theorem rests on the use of Boltzmann-Gibbs (BG) entropy. Therefore, it is useful to analyze transition 
regimes where the stationary distributions are of exponential form. On the other hand, the study of open systems 
may result in stationary distributions of the inverse power law form. Motivated by this fact, a generalization of 
the S-theorem has recently been provided by the employment of nonadditive Tsallis entropy [24, 25]. This gener- 
alized nonadditive S-theorem can be used to analyze the transition regimes in open systems, where the stationary 
distributions are g-exponentials, which asymptotically yield to inverse power law distributions. 
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Although the dissipative optical lattices serve as open systems due to their interaction with laser fields, neither 
original S-theorem due to Klimontovich nor its nonadditive generalization can be used for the description of self- 
organization in the aforementioned intermediate region of atomic transport. The reason is that this region exhibits 
a transition from Gaussian to g-Gaussian distributions. Therefore, the intermediate region requires an intermediate 
treatment between the two formulations of the S-theorem. This intermediate S-theorem, called modified S-theorem, 
is provided here for the explicit purpose of applying it to the dissipative optical lattices. 

The paper is organized as follows: In Section II, we review the original (i.e., additive) and nonadditive S-theorems 
for open systems. The modified S-theorem and its immediate application to dissipative open lattices are presented in 
Section III. Concluding remarks will be presented in Section IV. 



II. THE ADDITIVE AND NONADDITIVE S-THEOREMS 



The additive S-theorem due to Klimontovich [12-17] requires the use of two distinct probability distributions, i.e., r eq 
and p, corresponding to equilibrium and noncquilibrium states, respectively. The stationary equilibrium distribution 
is defined as the distribution corresponding to the state where the relevant control parameter is set to zero. Similarly, 
any other stationary state with non-zero control parameter is defined as the nonequilibrium state. As the value 
of the control parameter increases, the system recedes away from equilibrium state. By Gibbs' theorem, we know 
that mere substitution of these distributions into BG entropy will yield nonsensical results due to the failure of the 
equal mean energy assumption because of the interaction with the environment. Therefore, Klimontovich [12-17] 
renormalizes the equilibium distribution i.e., r eq . However, we note that the S-theorem of Klimontovich is not limited 
to the comparisons of equilibrium and nonequilibrium distributions. It can be used to compare any two stationary 
distributions as long as they are different from one another through the change in the value of the control parameter. 
From now on, we will drop the subscript (eq) from the equilibrium probability distribution r. Therefore, it should 
be understood that the probability distribution r denotes the equilibrium distribution. Next, Klimontovich defines a 
quantity called effective energy U e e [12-17] in terms of the equilibrium state as 



C/ eff =-lnr. (1) 

The definition of effective energy in Eq. (1) is central to the S-theorem, and therefore requires some explanation. The 
application of the effective energy definition in Eq. (1) to the equilibrium distribution results in U e e = (apart 
from normalization). This explains why it is called effective energy since it is proportional to the multiplication 
of the Lagrange multiplier (3 associated with the internal energy constraint and the energy of the zth microstate. 
Consequently, its average corresponds to the (3 times the average internal energy of the equilibrium distribution. We 
then equate the effective energy of the equilibrium and nonequilibrium states and this equalization procedure is called 
renormalization by Klimontovich. As a result of this renormalization, we no longer use the equilibrium distribution 
r, but the renormalized equilibrium distribution f, where tilde denotes that the renormalization is taken care of. 
Explicitly written, the renormalization reads 

(U cS ) ireq) = (t/ off > (ne9) , (2) 

where superscripts (req) and (neq) denote the renormalized equilibrium and ordinary noncquilibrium states, respec- 
tively. Then, S-theorem is equivalent to showing that the renormalized entropy i£[p||r] defined as 

R\p\\7\ = S^S(P) S%%(7) (3) 
is negative i.e., J?[p||r] < 0, since this implies that Sg G > S^q, where Sbg is the usual BG entropy given by 

w 

Sbg(p) = -^Pilript, (4) 

i=l 

where pi is the probability of the system in the ith microstate, W is the total number of the configurations of the 
system. 

The proof of the S-theorem relies on the definition of Kullback-Lcibler relative entropy (KL) [26], which reads 
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W , V 

K[p||r] = ^ Pi ln (f , (5) 

»=i 

since the renormalized entropy given in Eq. (3) is, through the renormalization in Eq. (1), equal to 

R\p\\r\ = -K\p\\r\. (6) 
The S-theorem is therefore proven [27], since KL entropy is positive definite i.e., 

R\p\\r\ < 0- (7) 

Although Klimontovich's S-theorem is a generalization of the Gibbs' theorem, it is based on the use of exponential 
stationary distributions and therefore BG entropy. Therefore, the S-thcorem due to Klimontovich cannot be used 
in the case of inverse power law stationary distributions. In order to overcome this difficulty, a new generalization 
of S-theorem has been given in the framework of nonadditive thermostatistics. This nonadditive S-theorem first 
generalizes the effective energy by deforming the logarithm in Eq. (1) as 

U cS = ln q (—) , (8) 
where the g-logarithm function ln g (x) is defined as 

r l-9 _ 1 

\n q {x) = (9) 
The nonadditive S-theorem can be similarly written as 

R q \p\\7\ = S q ( P ) S q (r) = -K q [p\m < 0, (10) 
for positive values of the nonadditivity parameter q, where the nonadditive Tsallis entropy S q [24, 25] is given by 

ir^W q _ i 

Sq(p) - ^\ f~ q (11) 



and the nonadditive relative entropy if g [p||r] [28] is given by 



K q 



The proof of the nonadditive S-theorem is limited to the positive q values only, since the nonadditive relative entropy 
is positive only in this region. However, we can exclude negative q values, since Tsallis entropy is stable only for 
positive values of q [29] . 

It should be noted that the ordinary S-theorem by Klimontovich is recovered in the q — > 1 limit. This can be seen 
from the inspection of Eq. (10), since Tsallis entropy expressions become BG entropies, whereas all the stationary 
distributions become exponential in the q — > 1 limit. The nonadditive relative entropy in Eq. (12) becomes KL 
relative entropy in the aforementioned limit i.e., if[p||r] = J2iPi^ n (Pi/ r i) [26, 28], which is positive definite, ensuring 
the negativity of the ordinary renormalized entropy. 
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III. OPTICAL LATTICES AND MODIFIED S-THEOREM 

In the previous section, we reviewed two interrelated criteria of self-organization, namely, the additive and non- 
additive S-theorems. The former can be used when there are exponential stationary distributions so that one can 
use BG entropy together with the renormalization of effective energies of stationary states. The latter, on the other 
hand, is used when the stationary distributions are g-exponentials associated with the Tsallis entropy. However, 
self-organization is not limited to the transitions amongst exponential or q-exponcntial stationary distributions. An 
intermediate transition is possible so that a transition from an exponential to ^-exponential stationary distribution 
may be observed. In this case, neither additive nor nonadditivc S-thcorems can be used. 

Such a transition has recently been observed in anomalous transport experiments in optical lattices [10]. Before 
proceeding further, it is important to review how the transition from an ordinary exponential stationary distribution 
to the q-exponcntial distribution can occur [9]. In order to do so, one first needs to derive the master equation, which 
describes the atom-laser interaction in the optical lattice. Assuming that the laser intensity is low and the atoms 
move fast, a semiclassical treatment is possible, which allows one to obtain, after spatial averaging, so-called Rayleigh 
equation for the semiclassical Wigner function W(p, t) 



dW d d 
_ = - [*(p)W] + s 



'nf \ dW 



(13) 



where the variables p and t denote momentum and time, respectively [30, 31, 32]. The drift coefficient K{p) and 
diffusion factor D(p) satisfy the relation [33] 

K{ P ) = (5 dU(p) 

D(p) ' l-[3(l-q)U(p) dp ' 1 ! 

where 

q=l + a c ; U (p) = p 2 . (15) 

The parameter j3 is positive and can be expressed in terms of the microscopic damping and diffusion coefficients [9] . 
The term a c denotes the control parameter, which represents the interaction of the atom with the environment i.e., 
laser field. The control parameter can be written more explicitly in terms of the microscopic parameters [9] as 

a c = — ( 16 ) 
api 

The terms a and p c denote the damping coefficient and the capture momentum, respectively. The term Dq is the 
diffusion term associated with the fluctuations due to spontaneous photon emissions and fluctuations in the difference 
of photons absorbed in the two laser beams [9]. Note that the control parameter a c can also be written in terms of 
the potential depth Uo and the recoil energy En [30] as 

a c=~F^- (17) 

Since the condition (14) is satisfied in optical lattices, the exact stationary solution to the Rayleigh equation in Eq. 
(13) is given in terms of the g-Gaussians 

W g (p) = C q [l-f3(l-q)U (p)] 1/{1 - q) , (18) 

where we also require W q (p) — ► when p — > ±oo When the control parameter is equal to zero, the nonadditivity 
parameter q becomes equal to 1. We then have stationary equilibrium distribution of the Gaussian form i.e., 



W 1 (p) = C 1 exp[-pU(p)} 



(19) 
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As the interaction between the atoms and the laser field is turned on, the control parameter a c takes nonzero positive 
values. As a result, the value of the nonadditivity parameter q exceeds 1 so that we have stationary nonequilibrium 
distribution of g-Gaussian form given by Eq. (18). 

As we remarked earlier, this transition from an ordinary exponential stationary distribution to a g-Gaussian dis- 
tribution cannot be analyzed within the frameworks of the aforementioned versions of the S-theorem. Therefore, we 
need a modified S-theorem, which can be used as a criterion of self-organization for the transition from Gaussian to 
g-Gaussian regime. 

Since the equilibrium distribution is an ordinary Gaussian, the modified S-theorem must include BG entropy asso- 
ciated with the Gaussian equilibrium distribution as is the case with the original, additive S-thcorcm by Klimontovich 
[12-17]. On the other hand, the stationary nonequilibrium distribution is a q-Gaussian. Therefore, the entropy of the 
stationary nonequilibrium distribution must be calculated by using Tsallis entropy as in the nonadditive S-theorem 
[24, 25]. The modified S-theorem will be the difference of these two entropies together with the renormalization 
condition. 

The only ingredient left is how the renormalization will be carried out in the modified S-theorem. However, 
the renormalization condition can be borrowed from Klimontovich's S-theorem [12-17] as it is, since the equilibrium 
distribution associated with the original and modified S-theorems both is of the same form i.e., exponential. Therefore, 
similar to Eq. (1), we can write 



U eS =-lnW 1 (p) 



(20) 



and then use Eq. (2) to obtain 



2\(reg) _ , 2\{ ne l) 



( P 2 r q ' = (p 2 ) 



(21) 



where we have inserted Eq. (20). The renormalized entropy SftfW^H W\] corresponding to the modified S-theorem can 
then be written as 



K[W q \\W 1 } = S q {W q )-S BG {W 1 ). 



(22) 



Similar to the additive and nonadditive versions of S-theorem, we calculate 3f[W g || W{\ and see whether it becomes 
negative for some interval of q values. Before proceeding with the renormalization condition given by Eq. (21), we 
write the equilibrium 



Wi(p) 

and the stationary nonequilibrium distribution 



- exp [-/3p 2 ] 



(23) 



W q (p) = y/4P(q-l) 



B 



1 2>-q 

2' 2q-2 



[l~(3{l-q)p 2 ] 



21 1/(1-9) 



; K«<3, 



(24) 



by inserting the normalization constants C\ and C q explicitly. The term B(x,y) denotes the Beta function [34]. We 
now rewrite Eq. (21) 



dpp 2 Wi (p) = J dpp 2 W q (p) . 



By substituting Eqs. (23) and (24) into Eq. (25), we obtain the renormalized inverse temperature (3 as 



(25) 



~ (3- 1 ) r (l^) fl , 5 



(26) 



G 



where T(x,y) is the Gamma function [34]. The associated equilibrium 



; ±<q< 



and nonequilibrium entropies can be calculated 



(27) 



s q to 



1 - [40 (q 1)] V [B (I, £=*)]"' B (§, J^) 5 
— ^ ; K?<g, 



9-1 



(28) 



where the renormalized inverse temperature in Eq. (27) can be substituted from Eq. (26). Finally, the rcnormalized 
entropy reads 



H[W,||Wi 



l-[Ap( q -l)]^[ B [l^)Y q B[\^) 
5-1 

• 4 ( g -i)r(^ 



1 i . 

- + - In 

2 2 



; K?< 



(29) 



The renormalized entropy 5?[W q ||Wi] in Eq. (29) is indeed negative in the interval 1 < q < | as can be seen from 
Fig. 1, independent of the value of 0. It corresponds to an entropy decrease when the entropies of the states are 
evaluated with the consistent expressions of entropy i.e., the Gaussian distribution at equilibrium with BG entropy 
and (/-Gaussian distributions at nonequilibrium with the Tsallis entropy. As the control parameter a c increases 
from zero (i.e., as the nonadditivity parameter q increases from 1), the system moves from Gaussian distribution to 
g-Gaussian distributions, decreasing its entropy. This indicates that the transition from the Gaussian distribution 
to the (/-Gaussian distributions in dissipative optical lattices is self-organization, signifying a transition to a more 
ordered state. It can be seen from Fig. 2b of Douglas et al. [10] that the self-organization interval determined by the 
modified S-theorem i.e., 1 < q < |, is confirmed by the experiment, although the whole range of admissible q values 
is theoretically 1 < q < 3. 

We also remark that the renormalized inverse temperature is a linear function of the ordinary inverse temper- 
ature as can be seen from Eq.(26). Therefore, the procedure of renormalization does not deform the equilibrium 
distribution, which is an exponential. Its role is to compensate the difference between the equilibrium temperature 
associated with the Gaussian distribution and the parameter in the (/-Gaussian distribution, since the latter is not 
physical temperature. This compensation is equivalent to heating up the equilibrium temperature of the system as a 
result of the linear decrease of with respect to for all q values in the range 1 < q < |. The fact that depends 
explicitly on q implies that the renormalization also takes the microscopic structure into account. 

We finally mention the effect of renormalization on the spatial correlation function G (a) = J dxtp (x, t) ip* (x + a,t), 
where ip (x,t) denotes the wave function at time t [35]. The spatial correlation function can be considered as a measure 
of how coherent a state is between two different points. For a Gaussian wave packet such as W\, the spatial correlation 
function too is a Gaussian with a correlation length A = 2hy/0 [9]. The renormalized correlation length, on the other 



hand, is equal to A = 2hJ (q — l)T (j^j /T in the interval 1 < q 



< 



Therefore, we conclude that the 
In other words, the Gaussian 



2q-2 J I \2q 

correlation length decreases with increasing value of the nonadditivity parameter q. 
distribution becomes less coherent spatially as a result of the renormalization procedure employed in the modified 
S-theorem. This decrease in coherence becomes more dominant as the system moves further from equilibrium. 



IV. CONCLUSIONS 



By analyzing the intermediate anomalous transport regime in dissipative optical lattices [9, 10] in the framework 
of a modified version of Klimontovich's S-theorem, we have shown that the self-organization occurs in the transition 
from Gaussian to the (/-Gaussian distributions, implying an entropy decrease i.e., a more ordered state as the system 
is driven out of equilibrium as a result of the interaction of the atoms with the laser field. 
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An interesting result of the modified S-theorem is that it confines the self-organization in the region 1 < q < |, 
although, for example, there is nothing which forbids transition to a region with q = 2, (note that the q-Gaussian 
distributions are normalizable in the whole range 1 < q < 3). This is also confirmed with the experimental results 
obtained by Douglas et al. [10]. The reason for this confinement is due to the additional requirement associated with 
the (modified) S-theorem i.e., renormalization of the effective energy. For dissipative optical lattices, the effective 
energy term corresponds to the square of the momentum. Therefore, the procedure of renormalization confines the 
self-organization to the region, where the second moment of the distribution (p 2 ) converges. This is similar to the 
requirement of finite capture momentum for the occurrence of g-Gaussian distribution in dissipative optical lattices 
[9]. If the capture momentum p c is infinite, then the control parameter a c becomes zero according to Eq. (16). Since 
zero value of the control parameter corresponds to the case q = 1, we obtain only Gaussian dynamics, eliminating the 
possibility of an anomalous transition [9]. The existence of the non-Gaussian transition requires both finite capture 
momentum p c and finite second moment (p 2 )- 

We also conclude that the correlation length associated with the equilibrium distribution decreases with increasing 
value of the nonadditivity parameter q as a result of the renormalization procedure. This implies that the equilibrium 
Gaussian distribution becomes less coherent spatially as the system moves further from equilibrium. 

We finally note that another version of the modified S-theorem was provided by Bashkirov and Vityazev in Ref. 
[36] based on Rcnyi entropy [37]. However, the compared distributions were both chosen as the inverse power law 
distributions although the renormalized entropy in Ref. [36] was defined as the difference of Renyi and BG entropies. 
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FIG. 1: The renormalized entropy SR[W q ||VVi] in Eq. (29) is plotted in the interval 1 < q < | for a representative range of (5, 
from which the negativity of the renormalized entropy is evident. 



